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JJ^ Here we study the static potential for the two quarks and two antiquarks system. First this is done 

using the tetraquark operator, which has been previously used to calculate the static potential. 
This is found to give good results in the region where the tetraquark is expected to be the ground 
state, however failing outside it. To repair this, we resort to a variational method. This let us study 
the first excited state besides the ground state of the system in two different particle dispositions, 
one where the quarks are on the same side of a rectangle and the other where they are at opposite 
sides. Results for the field components and for the lagrangian density are presented. 
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Figure 1: Wilson loop operator used by [2] and [3] to calculate the static potential of a tetraquark system. 

1. Motivation 

Systems constituted by two quarks and two antiquarks are of extreme importance for strong 
interaction physics. Not only because they are a starting point for meson-meson scattering, but 
also because of the possible existence of bound-states — tetraquarks, initially predicted by Jaffe 
[1]. There are several observed resonances who are candidates to tetraquarks. The most recent 
candidates are the particles reported by the Belle collaboration. Recently, static systems of two 
quarks and two antiquarks, were studied in the lattice [2, 3] using an operator with the quantum 
numbers of the tetraquark. The results indicate that the static potential of this system could be well 
described by the generalized flip-flop potential. This potential, is similar to the ones which were 
first proposed [4, 5] as a device to suppress the non-physical Van der Waals forces, which arise 
when potential models based on Casimir Scaling are used, differing by the existence of a third 
branch of the potential — the tetraquark branch. 

2. Tetraquark Wilson Loop 

The Static potential for the tetraquark has been studied in the lattice by [2] and [3]. 
For that, a Wilson Loop operator was constructed. It is given by 

W4 Q = ^Tr[M 1 R 12 M 2 L 12 i (2.1) 

with 

Rn = ei jk £i>fk>R{ f R k 2 k ' ( 2 - 2 ) 

j ii' e c T jj' T kk! 

^12 — ^ijk^i' j'k ,l ^\ ^2 ■ 

The corresponding paths are visually described on figure 2. 

This Wilson Loop has the quantum numbers of colour singlet system where the two quarks 
form an antitriplet and the two antiquarks form a triplet. 



2 



Tetraquark and the flux tube recombination 



Marco Cardoso 



A 




Figure 2: Here we can see the shape of the minimal string which links the four particles in the tetraquark. 
Note that when the diquark and the diantiquark are close the five segment structure collapses into a four 
segment one. 



This lattice studies, indicate that the static potential of the two-quark and two-antiquark system 
is a generalized flip-flop potential: 

Vff min(Vr, V Mi m 2 ,Vm 3 m 4 ) (2.3) 

where Vm x m 2 an d ^m 3 m 4 are the two possible two-meson potentials, given by the sum of two in- 
dependent intra-meson potentials VmiM 2 = Vm 1 +Mm 2 , and Vj is the tetraquark potential which 
corresponds to the sector where the four particles are confined, linked by a single fundamental 
string. It is given by 

X- X ■ 

V r = C+ a£-^ + oL min , (2.4) 

i<j 1 1 

where L m in is the minimal distance linking the four particles, which corresponds to the string 
configuration displayed on fig. 2. 

For the special case where the four particles form as rectangle as in fig. 2, L min is given by 
Lmin = V^n + r 2 for r 2 > ^ [6, 7]. If we neglect the Coulomb part of Vff, we could estimate that 
the tetraquark sector becomes the ground state when r 2 > V^n . 



3. Tetraquark Fields 

The fields in the tetraquark system can be calculated by using the correlation of plaquettes and 
Wilson Loops. The Squared chromoelectric and chromomagnetic fields, are computed by 

<*?> = <>>»,> (3-D 
where the index i complments j and k. P^ v is P^ v — I — jTr[[/ jU (s)[/ v (s + v)[/^(s + v)ul(s)]. 
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Figure 3: Plot of the lagrangian density for different values of r\ with r2 = 14. 
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Figure 4: Comparision between the squared chromoelectric field in the quark- antiquark flux tube core (in 
the meson) and the diquark-antidiquark flux tube core (in the tetraquark). 



The Lagrangian and energy densities are given by 

(if) = \(E 2 -B 2 ) (3.2) 
<JT) = \{E 2 + B 2 ). 

Now, we will present results for the chromo-fields for the tetraquark system for the simple 
case where the four particles form a rectangles, as in fig. 2. 

This results, and the one of the following chapters were obtained using quenched lattice QCD 
configurations with dimension 24 3 x 48 and j8 = 6.2. They where used in GPUs by using a com- 
bination of Cabbibo-Marinari, pseudo heat-bath and over-relaxation algorithms [8]. APE smearing 
[9] and Hypercubic blocking [10] were used to improve the signal to noise ratio. 

The results for the lagrangian density are given on fig. 3 for fixed r2 = 14. 

As we can see the in the leftmost fig. r\ = 0, in which case the system the system is collapsed 
in to a meson. Then we can see left to right, with the increase of r\, the confining string acquires a 
shape similar with the one on fig. 2 for r2 > in agreement with the result obtained for the static 
potential. 
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Figure 5: Lagrangian density in the flux tube for r\ = r2 = 8 and r\ = 8, r2 = 14. 



In fig. 3 it is compared the squared chromoelectric fields (E 2 ) in the flux tube centre of the 
meson and in the centre of the diquark-diantiquark flux tube in the tetraquark. As can be seen 
both the flux-tubes have a similar behaviour, which confirms that the string which is present in this 
system is a fundamental one, again in agreement with the results for the static potential. 

However, as can be seen in fig. 3 this agreement disappears for Y2 < \fi>r\ — the region where 
the ground state of the system is a two meson state and not the tetraquark state. As can be seen there, 
for r\ = 8 and ri = 14 (tetraquark region), we obtained the expected ground state. However for 
r\ = r2 = 8 (meson-meson region) we still see a string linking the four particles as in the tetraquark 
state, but we should see, to be in agreement with the generalized flip-flop model supported by the 
results obtained for the static potential, two strings corresponding to the two mesons. We think 
that this disagreement is due to the low overlap of the used operator with the ground state, and also 
to the small temporal extensions used in the Wilson Loop. To obtain the true ground state of the 
system in this region and also the first excited state we used a variational method. 

4. Variational Method for the QQQQ system 

Now we will improve this method in order to obtain the true ground state of the system and 
also to obtain the first excited state. To achieve this, we note that the Wilson Loop operator could 
be written as a correlation of a certain operator at different times W(t) = {@(t)ffi(0)). This could 
be generalized by considering instead a base of operators &i. So this way, our wilson loop becomes 
a matrix = (^^J). This could be used not only to improve the ground state overlap but also to 
obtain more energy levels of the system. These are given by solving the generalized eigensystem 

(W ij (t))4(t)=w n (W ij (0))4(t) (4.1) 

The fields can then be calculated by (P) n = — (P) with W n = c l n WijC J n . 

We will consider again the case where the four particles form a rectangles and also two differ- 
ent alignments of the QQQQ system. A parallel one, where the two quarks are on the same side 
of the rectangle and an antiparallel one, where the quarks are on opposite corners of the rectangle 
(fig. 4). 

For both cases we will use a base of two operators, similarly to what was made in [11]. In 
the parallel case the operators will be the tetraquark operator &$q, the correlation of which is the 
tetraquark Wilson Loop W^q = (^4g(0^4g(0))» an( i a two meson operator. This gives a Wilson 
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Figure 6: Left: Parallel alignment. Right: Antiparallel alignment. 
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Figure 7: Elements of the Wilson loop matrix used for the parallel alignment. 



Figure 8: Elements of the Wilson loop matrix used for the antiparallel alignment. 

loop matrix where the diagonal elements are W4Q and the correlation of two wilson loops, while 
the off-diagonal elements correspond to the transition between the two states (see fig. 4 

For the antiparallel alignment, the operators used were the two meson-meson operators, giving 
the four matrix elements given on fig. 4. 

5. Results 



Now, we present the results obtained for the parallel geometry. 

In fig. 5 we can see the components of the chromoelectric and chromomagnetic fields for the 
ground state of the parallel alignment, with r\ = 6 and ri = 12, while in fig. 5 the componentes are 
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Figure 9: Components of the chromo-electromagnetic field, for the ground state of the parallel geometry 
with r\ = 6 and 7*2 = 12. 



presented for the ground state of the antiparallel alignment with n = 6 and 7*2 = 10. The results 
indicate that what we have is essencialy a tetraquark in the first case and a two meson system in the 
second. Besides that we can see features that are common to flux-tubes in general [9] such as the 
dominance of the longitudinal chromoelectric field. 

In fig. 5 it is shown the evolution of the lagrangian density of the ground state of the parallel 
geometry, for fixed r\ = 6. Here we can see the transition between the two meson state in the top- 
left picture (7*2 = 6) to the tetraquark state in the bottom right picture 7*2 = 12. This is in agreement 
with the predictions of the generalized flip-flop potential contrary to the results we obtained using 
only the tetraquark operator W4Q. The same evolution is presented in 5, but for the first excited 
state of the same configurations. This results are not so readily understandable. 

To better compare the ground and the first excited state we can look at fig. 5 in the y = axis. 

In fig. 5 we see the lagrangian density for the parallel alignment with varying 7*1 and 7-2- For 
7*1 = 6 and 7*2 = 10 we can see the two mesons almost completely separatedly, while for r\ = 7*2 
we have a state when all the four particles are linked. We only show pictures for 7*2 > n . The case 
7*2 < 7*1 can be seen by rotating the pictures for 7*2 > r\ . This behaviour is again expected, as we 
obtain the separation in two mesons of the system when it correspond to the minimal string energy. 
The behaviour for r\ = 7*2 is also expected, because neither of the two meson systems is preferred 
over the other. 

The results for the first excited state are presented in fig. 5. As we can see there, the first excited 
state is different from a two meson state, with all the four particles being linked the flux-tube. 

Finally we compare the ground state and the first excited state in this geometry in fig. 5 for 
y = 0, in the left, for r\ = 6 and 7*2 = 10, in the right for 7*1 = 7*2 = 10. As can be seen, for the 
first case, the ground state flux-tube in this region is very small, because of the formation of two 
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Figure 10: Components of the chromo-electromagnetic field, for the ground state of the anti-parallel geom- 
etry with r\ = 6 and Y2 = 10. 
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Figure 11: Lagrangian density for different values of r\ and r2 in the ground state of the parallel geometry. 
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Figure 13: Comparasion of the action density in the center of flux tube in the parallel geometry for r\ = 6 
and Y2 — 12. 



mesons. 



6. Discussion 



Our results for the fields of the two quarks and two antiquarks system agree with the results 
obtained for the static potential, therefore supporting a string picture of confinement with flux-tube 
recombination. Besides this, we were also able to calculate the first excited state of this system. 
The low overlap between the tetraquark operator and the two meson ground state may indicate 
small width tetraquark resonances. We are also studying the fields of the pentaquarks and the static 
potentials of this system. 
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Figure 14: Lagrangian density for different values of r\ and in the ground state of the anti-parallel 
geometry. 
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Figure 15: Lagrangian density for different values of r\ and in the first excited state of the anti-parallel 
geometry. 
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Figure 16: Comparasion of the action density in the center of flux tube in the anti-parallel geometry for two 
geometries. 
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